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ABSTRACT 

The deep, wide-area (~800-900 arcmin 2 ) near-infrared/WFC3/IR + Spitzer /IRAC observations over 
the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift 
UV luminosity functions (LFs). However, the lack of HST 1.05/jm observations over the CANDELS 
fields has made it difficult to identify z ~ 9-10 sources robustly, since such data are needed to confirm 
the presence of an abrupt Lyman break at 1.2 /irn. We report here on the successful identification 
of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest- 
probability z~9-10 galaxy candidates with observations at 1.05/xm, to search for a robust Lyman- 
break at 1.2p.m. The potential z~9-10 candidates are preselected from the full HST, Spitzer /IRAC 
S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations 
(CFHTLS-DEEP+HUGS+UltraVISTA+ZFOURGE). We identified 15 credible z~9-10 galaxies over 
the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted 
follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a 
reliable sample of z ~ 9-10 galaxies. Through extensive simulations, we replicate the selection process 
for our sample (both the preselection and follow-up) and use it to improve current estimates for 
the volume density of bright z^9 and z~10 galaxies. The volume densities we find are 5+1 x and 
8l|x lower, respectively, than found at z ~ 8. When compared with the best-fit evolution (i.e., 
dlog 10 puv/dz = —0.29 ± 0.02) in the UV luminosities densities from z ~ 8 to z ~ 4 integrated to 
0.3L*_3 (—20 mag), these luminosity densities are 2.6 +q |)x and 2.2+1°x lower, respectively, than the 
extrapolated trends. Our new results are broadly consistent with the “accelerated evolution” scenario 
at z > 8, as seen in many theoretical models. 

1. INTRODUCTION 

The first galaxies are thought to have formed in the 
first 300-400 Myr of the universe. Over the last few years, 
remarkable progress has been made in extending samples 
back to this time, with more than ~700 probable galaxies 
identified at z > 6.3 with HST (Bouwens et al. 2015; 

McLure et al. 2013; Finkelstein et al. 2015) and 20-30 
candidate galaxies identified as far back as redshifts 2 ~ 

9-11 (Bouwens et al. 2011a, 2014b; Zheng et al. 2012; 

Coe et al. 2013; Ellis et al. 2013; McLure et al. 2013; 

Oesch et al. 2013, 2014, 2015; Zitrin et al. 2014; Ishigaki 
et al. 2015; Bouwens et al. 2015; McLeod et al. 2015). 
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At present and over the next year, considerable re¬ 
sources are being devoted to both the discovery and 
study of ultra-faint galaxies with HST from the new 
Frontier Fields initiative (e.g., Lotz et al. 2014; Coe et 
al. 20150 The goal of this initiative is to combine the 
power of gravitational lensing from galaxy clusters with 
very deep exposures with the Hubble and Spitzer Space 
Telescopes. 840 orbits of HST observations are being 
invested in deep optical/ACS + near-IR/WFC3/IR ob¬ 
servations of 6 galaxy clusters. Deep observations of a 
“blank” field outside the galaxy clusters are also being 
obtained in parallel with observations over the clusters. 

Despite the considerable focus by the community on 
the Hubble Frontier Fields observations over galaxy clus¬ 
ters and deep fields (e.g., Atek et al. 2014, 2015; Zheng et 
al. 2014; Coe et al. 2015; Ishigaki et al. 2015; Oesch et al. 
2015), it is also possible to uncover modest numbers of lu¬ 
minous 2 ~ 9-10 galaxies over wide-held surveys as Oesch 
et al. (2014) first illustrated through the identification 
of six intrinsically-luminous 2 ~ 9-10 candidate galaxies 
over the GOODS-North and GOODS-South CANDELS 
fields (Grogin et al. 2011; Koekemoer et al. 2011). These 
sources allowed us to set some initial constraints on the 
rate at which t/U-luminous galaxies evolve with cosmic 
time, as well as providing some constraints on the ap¬ 
proximate shape of the UV LF at z = 9-10. 

Due to the inherent brightness of such sources, these 

10 http://www.stsci.edu/hst/campaigns/frontier-fields/ 


2 


sources are also valuable for efforts to measure the phys¬ 
ical properties of galaxies at very early times. Measure¬ 
ments of the DU-continuum slopes (Oesch et al. 2014; 
Wilkins et al. 2016), Balmer-break amplitudes (Oesch et 
al. 2014), stellar masses (Oesch et al. 2014), and sizes 
(Holwerda et al. 2014; Shibuya et al. 2015) can all be 
achieved using very bright galaxies. 

Despite the usefulness of bright 2 ~ 9-10 galaxies for 
addressing many contemporary science questions, cur¬ 
rent samples of these objects remain quite small and 
constraints on their volume densities poor. Of partic¬ 
ular note, the recent Oesch et al. (2014) sample over the 
GOODS-North and GOODS-South fields only contained 
2 bright z ~ 9 and 4 bright z ~ 10 galaxies. With such 
small samples, current uncertainties on the volume den¬ 
sity of bright 2 ~ 9-10 galaxies are large indeed (> 0.3- 
0.4 dex). This is especially the case when one considers 
the impact of field-to-field variations (“cosmic variance”) 
which is as large as a factor of two across the CANDELS 
fields, e.g., see Figure 14 from Bouwens et al. (2015), and 
may be even larger for the brightest sources (Bowler et al. 
2015; Roberts-Borsani et al. 2016). Clearly, we require 
many independent lines of sight on the z ~ 9-10 universe 
to average over the large-scale structure. Unfortunately, 
the Frontier Fields Initiative will not significantly help 
with this issue for luminous sources, given the limited 
area covered by observations from this program. 

Nevertheless, there is a huge quantity of HST and 
Spitzer data already available that can be used to con¬ 
struct larger samples of bright z ~ 9-10 galaxies. The 
most significant of these data sets are the ~500 arcmin 2 
CANDELS UDS, COSMOS, and EGS fields which fea¬ 
ture very deep optical, near-IR, and Spitzer /IRAC obser¬ 
vations. These observations are very useful for the robust 
detection of bright z ~ 9-10 candidates and also confirm¬ 
ing a blue color redward of the break, distinguishing such 
galaxies from dusty, red galaxies at z ~ 1-3. While pos¬ 
sessing great potential, the CANDELS UDS, COSMOS, 
and EGS search areas lack correspondingly deep obser¬ 
vations at 1.05/um, just blueward of the Lyman-break in 
candidate z ~ 9-10 galaxies which is important for con¬ 
firming a spectral break at ~1.2/zm and also distinguish¬ 
ing these z > 9 galaxy candidates from Balmer-break 
sources at 2 ~ 1-3. 

Fortunately, we can overcome the aforementioned lim¬ 
itations of the CANDELS UDS, COSMOS, and EGS 
data sets by leveraging essentially all existing observa¬ 
tions over these fields (Table U) to first identify the high¬ 
est probability z ~ 9-10 candidates over these fields and 
then obtaining targeted follow-up observations of these 
candidates at 1.05/zm to determine which are likely at 
2 > 8 (Figure [TJ. In cycle 21, we successfully pro¬ 
posed for such a follow-up program of plausible candidate 
2 ~ 9-10 galaxies over the CANDELS UDS, COSMOS, 
and EGS fields. Observations from this program - which 
we call z9 (redshift 9)-CANDELS (Bouwens 2014: GO 
13792) are now complete and cover all 12 of the primary 
candidates from that program. Based on the informa¬ 
tion we obtain from our proposed follow-up observations 
and the selection criteria we use in identifying our initial 
sample of 12 candidate 2 ~ 9-10 galaxies from these three 
CANDELS fields, we can derive constrants on the volume 
density of luminous 2 ~ 9-10 galaxies. Searches over 
CANDELS-GOODS-North, CANDELS-GOODS-South, 


the ERS fields can be further used to improve the con¬ 
straints we can obtain on the bright end of the z ~ 9-10 
LFs. 

In this paper, we describe (1) the preselection we uti¬ 
lized to identify candidate 2 ~ 9-10 galaxies from the 
CANDELS-UDS, COSMOS, and EGS fields for HST 
follow-up observations and (2) the results from this pro¬ 
gram. Our primary scientific objective is in obtaining the 
best available constraints on the volume density of espe¬ 
cially luminous z ~ 9 and 2 ~ 10 galaxies. Through such 
constraints, we have a direct measure of how fast (1) the 
bright end of the UV LF and (2) [TV-luminous galaxies 
evolve. Through comparison with the volume density of 
fainter sources, the present search results also allow us 
to constrain the overall shape of the UV LF. Finally, we 
would expect our selection to allow us to considerably 
expand the overall sample of bright 2 ~ 9 and 2 ~ 10 
galaxies available over the CANDELS fields. This has 
value both for the further characterization of the phys¬ 
ical properties of 2 > 9 galaxies and as bright sources 
to target with early JWST observations. These bright 
samples will be further enhanced with bright 2 ~ 9-10 
galaxies from the BoRG[ 2 gi 0 ] program (Trenti 2014). 

Here is a brief plan for this paper. In §2, we include 
a description of the observational data that we utilize 
to identify high-probability 2 ~ 9-10 galaxies over the 
CANDELS-UDS, COSMOS, and EGS fields. In §3, we 
describe our criteria for performing photometry and iden¬ 
tifying high-probability 2 ~ 9-10 galaxy candidates over 
the CANDELS-UDS, COSMOS, and EGS fields. In §4, 
we describe the results of the z9-CANDELS program 
where we use these observations to ascertain the likely 
nature of our selected 2 ~ 9-10 candidate galaxies. In §5, 
we describe our search results for bright 2 > 8.4 over the 
CANDELS GOODS-North, GOODS-South fields, and 
ERS fields, extending previous work by Oesch et al. 
(2014: see also McLure et al. 2013 who also conducted 
such a search over the GOODS-South field). Finally, in 
§6, we make use of these search results to provide the first 
constraints on the bright end of the 2 ~ 9 and 2 ~ 10 
LFs using a search over all five CANDELS fields. 

For consistency with previous work, we quote results 
in terms of the luminosity L* =3 Steidel et al. (1999) de¬ 
rived at 2 ~ 3, i.e., Mnoo ab = —21.07. We refer to the 
HST F435W, F606W, F775W, F814W, F850LP, F098M 
F105W, F125W, F140W, and F160W bands as R 435 , 

^606, *775) Is 14 , 2850 ) U)98) ^105> Jl25> JH 140, and -ffi60> 
respectively, for simplicity. Where necessary, we assume 
D 0 = 0.3, Da = 0.7, and H 0 = 70km/s/Mpc. All mag¬ 
nitudes are in the AB system (Oke & Gunn 1983). 

2. OBSERVATIONAL DATA 

In the present analysis, we conduct a search for bright 
2 ~ 9-10 candidate galaxies over the ~450 arcmin 2 re¬ 
gion within the CANDELS-UDS, COSMOS, and EGS 
fields with the deepest HST optical/ACS and near- 
IR/'WFC/IR observations (~75-80% of the WFC3/IR 
area). 

In conducting this search, we utilize the reductions of 
the HST observations described in Bouwens et al. (2015). 
Those reductions include all observations associated with 
the AEGIS, COSMOS, and CANDELS HST surveys 
and SNe follow-up programs, including the J-ffi 4 o-band 
observations associated with 3D -HST (Brammer et al. 
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TABLE 1 

Observational Data Used to Identify! the Bright z ~ 9-10 Candidate Galaxies over the CANDELS UDS, COSMOS, and EGS 

fields.* 


CANDELS UDS 

Two-Part Search Strategy (Preselection + Follow-up: 

CANDELS COSMOS 

§3, §4) 

CANDELS EGS 



5 (7 



5(7 



5(7 


Filter! 

Depth a 

Source 

Filter! 

Depth 3, 

Source 

Filter! 

Depth a 

Source 

+606 

26.8 

HST/ ACS 

+606 

26.5 

HST/ ACS 

+606 

27.3 

HST/ACS 

+814 

26.8 

HST/ ACS 

+814 

26.5 

HST/ ACS 

+814 

27.1 

HST/ ACS 

+1 25 

26.3 

++ST/WFC3 

+125 

26.1 

FST/WFC3 

+125 

26.4 

HST/ WFC3 

JH 140 

26.1 

++ST/WFC3 

JH 140 

25.8 

HST/ WFC3 

JH 140 

25.6 

++ST/WFC3 

++160 

26.5 

++ST/WFC3 

H 160 

26.3 

+LST/WFC3 

++160 

26.6 

HST/ WFC3 

u 

25.8 

CFHT/Megac.am 

u 

27.7 

CFHT/Megacam 

u 

27.4 

CFHT/Megacam 

B 

28.0 

Subaru / Suprime-Cam 

B + g 

28.4 

Subaru/Suprime-Cam + 

9 

27.8 

CFHT/Megacam 

V + r 

28.0 

Subaru / Suprime-Cam 



CFHT/Megacam 

r 

27.6 

CFHT/Megacam 

i 

27.4 

Subaru / Suprime-Cam 

V + r 

27.9 

Subaru/Suprime-Cam + 

i + y 

27.4 

CFHT/Megacam 

z 

26.3 

Subaru / Suprime-Cam 



CFHT/Megacam 

z 

26.0 

CFHT/Megacam 

Y 

25.9 

VLT/HAWKI/HUGS 

i + y 

27.7 

Subaru/Suprime-Cam + 

K 

24.1 

UKIRT / WIRCam 

+1 

25.6 

Magellan/FOURSTAR 



CFHT/Megacam 

3.6pm 

25.4 

Spitzer/S-CANDELS 

+2 

25.7 

Magellan/FOURSTAR 

z 

26.4 

Subaru/Suprime-Cam + 

4.5 pm 

25.3 

Spitzer/S-CANDELS 

J 

25.4 

UKIRT/WFC AM 



CFHT/Megacam 




+3 

25.4 

Magellan/FOURSTAR 

Y 

26.1 

UltraVISTA 




H 

24.6 

UKIRT/WFC AM 

Jl 

25.6 

Magellan/FOURSTAR 




H s 

25.0 

Magellan/FOURSTAR 

+2 

25.5 

Magellan/FOURSTAR 




++; 

24.8 

Magellan/FOURSTAR 

J 

25.3 

UltraVISTA 




K a 

25.5 

VLT/HAWKI/HUGS + 

+3 

25.3 

Magellan/FOURSTAR 






UKIRT/WFC AM + 

H a 

24.7 

Magellan/FOURSTAR 






Magellan/FOURSTAR 

H 

25.0 

UltraVISTA 




3.6pm 

25.4 

Spitzer/S-CANDELS 

Hi 

24.7 

Magellan/FOURSTAR 




4.5/.tm 

25.4 

Spitzer/S-CANDELS 

Ks 

25.3 

UltraVISTA + 









Magellan/FOURSTAR 







3.6pm 

25.3 

Spitzer/S-CANDELS 







4.5 pm 

25.3 

Spitzer/S-CANDELS 






Direct Search Strategy for z > 8.4 Galaxies (§5) 




CANDELS GOODS-South 

ERS 



CANDELS GOODS-North 

++435 

27.1-27.3 

HST/ ACS 

B435 

27.1 

HST/ACS 

B435 

27.2-27.3 

HST/ACS 

+606 

27.4-27.7 

HST/ ACS 

+606 

27.4 

HST/ ACS 

+606 

27.4 

HST/ ACS 

*775 + 



*775 + 



*775 + 



Is 14 

27.5-27.6 

HST/ ACS 

+814 

27.3 

HST/ACS 

+814 

27.2-27.7 

HST/ACS 

^850 

26.8-26.9 

HST/ ACS 

Z850 

26.7 

HST/ ACS 

^850 

26.9-27.0 

HST/ ACS 

1*105 

26.4-27.0 

++ST/WFC3 

1 098 

26.5 

HST/ WFC3 

+105 

26.5-26.8 

++ST/WFC3 

J 125 

26.5-27.0 

++ST/WFC3 

+125 

27.0 

HST/ WFC3 

+125 

26.4-27.2 

++ST/WFC3 

JH 140 

26.1 

++ST/WFC3 

JH 140 

25.8 

HST/WFC3 

JH 140 

25.6 

++ST/WFC3 

++ 160 

26.5-27.0 

HST/ WFC3 

Hieo 

26.9 

HST/ WFC3 

++160 

26.5-27.1 

++ST/WFC3 

K s 

26.5 

VLT/HAWKI/HUGS + 

Ks 

26.5 

VLT/HAWKI/HUGS + 






VLT/ISAAC + 



VLT/ISAAC + 






PANIC + 



PANIC + 






Magellan/FOURSTAR 



Magellan/FOURSTAR 




3.6pm 

25.8 

Spitzer/S-CANDELS 

3.6pm 

25.8 

Spitzer/S-CANDELS 

3.6pm 

25.8 

Spitzer/S-CANDELS 

4.5pm 

25.8 

Spitzer/S-CANDELS 

4.5 pm 

25.8 

Spitzer/S-CANDELS 

4.5 pm 

25.8 

Spitzer/S-CANDELS 


t For each source in our search fields, flux measurements are derived based on the all the observational data presented in this table. All of 
these measurements are used in deriving a redshift likelihood distribution for individual sources. 

a The 5 <j depth are estimated from the median 5(7 uncertainties on the total flux measurements of sources found over our search fields with 
/ 1 1 ,,(] ,\ g-band magnitudes of 26.0-26.5. 


2012) and AGHAST (Weiner et al. 2014) programs. 
Beyond the HST observations themselves, perhaps the 

most valuable data set that we can leverage in our 
search for probable z ~ 9-10 galaxies is the very deep 
Spitzer /IRAC S-CANDELS observations over the CAN- 
DELS fields (Ashby et al. 2015), which when combined 
with Spitzer/IRAC SEDS observations (Ashby et al. 

2013) reach 50 hours in depth (26.0 mag at 5<r: 2"- 
diameter apertures). Those observations provide us with 
constraints on the spectral slope of galaxies redward of 
the H i 6 o band, which when combined with evidence for 


a break across the J 125 and i+i 6 o bands and a non¬ 
detection at optical wavelengths is strongly suggestive 
of a z ~ 9-10 galaxy. 

In addition, we also make use of all significant, 
public ground-based observations over these fields, 
including optical observations from Subaru Suprime- 
Cam [CANDELS-COSMOS; CANDELS-UDS] and 
CFHT/Megacam [CANDELS-COSMOS; CANDELS- 
EGS] and deep near-IR observations from VISTA 
[CANDELS-COSMOS], UKIRT/WFCAM [CANDELS- 
UDS], and VLT/HAWKI [CANDELS-UDS], Magel- 
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tions over the COSMOS and EGS fieldsEB the public 
v2.0 reductions of the UltraVISTA observations (Mc¬ 
Cracken et al. 2012), the Cirasuolo et al. (2010) redutions 
of the deep Subaru Suprime-Cam observations over the 
UDS/SXDS field (Furusawa et al. 2008), the Bouwens 
et al. (2015) reductions of the HUGS HAWK-I obser¬ 
vations (Fontana et al. 2014), the public reductions of 
the WIRCam deep survey 7\ s -band observations over the 
CANDELS EGS field (McCracken et al. 2010; Bielby et 
al. 2012), the vO.9.3/vO.95.5 reductions of the ZFOURGE 
COSMOS/UDS observations (I. Labbe et al. 2015, in 
prep), the IUDF reductions of the Spitzer /IRAC observa¬ 
tions over the GOODS-South and GOODS-North fields 
(Labbe et al. 2015), and the public reductions of the 
Spitzer SEDS and S-CANDELS programs (Ashby et al. 
2013; Ashby et al. 2015). 

Table [T] provides a convenient summary of all the ob¬ 
servational data we utilize. Combining the flux measure¬ 
ments from the different data sets, the 5cr depths of these 
fields (derived from the median uncertainties on the to¬ 
tal flux measurements) range from ~28 mag at <0.8/rm, 
~26.0-26.5 mag at ~0.9/im, ~26.0 mag at 1.05/rm, ~26.6 
mag at ~1.2-1.6/tm, 24.1-25.5 mag at 2.3/rm, and 25.6- 
25.9 mag at 3.6^m + 4.5^m. We refer the interested 
reader to Figure 3 from Bouwens et al. (2015) for a graph¬ 
ical representation of these depths as a function of wave¬ 
length. 

3. Z ~ 9-10 SELECTION 


Wavelength (^m) 

Fig. 1.— (upper panel) Wavelength sensitivity curves for the 
filters (F606W, F814W, F125W, F160W) in which deep HST ob¬ 
servations are available over the CANDELS UDS, COSMOS, and 
EGS fields (black), as well as those (F105W) primarily obtained 
by our follow-up program z9-CANDELS (red), (lower panel) 5cr 
median depths of the observations versus wavelength available over 
the CANDELS UDS, COSMOS, and EGS fields (see also the com¬ 
pilations in Table 1 and Figure 2 of Bouwens et al. 2015). The 
depths plotted here are binned in such a way to combine all the 
data (Table [TJ that exist in 0.1-0.15jum segments. The depths do 
not include the ZFOURGE observations here, since those observa¬ 
tions only cover >65% of each CANDELS field (adding some useful 
depth from 1.0/xm to 1.7/zm). There is a modest wavelength gap 
between the deeper observations at ~ 0.6-0.9/xm and that which ex¬ 
ists at ~ 1.2-1.6/im. While some (~26-mag, 5 cr) observations exist 
at 1.05/xm to probe below the putative Lyman-break for z ~ 9-10 
galaxy candidates, the addition of deep observations at l^rni with 
HST can significantly improve current constraints on the robust¬ 
ness of the Lyman-break at lpm. 


lan/FOURSTAR [CANDELS-COSMOS; CANDELS- 
UDS], and CFHT/WIRCam [CANDELS-EGS], The 
deep optical observations allow us to search for faint 
optical flux in the z ~ 9-10 candidate identified over 
the CANDELS-UDS/COSMOS/EGS fields, while the 
near-IR observations allow us to test for the presence 
of a putative break at 1.2/im, verify that candidates 
show no flux blueward of the break, and test for a flat 
UV -continuum redward of the break. 

In our analysis of data over the five CANDELS fields 
and ~40 arcmin 2 ERS field (Windhorst et al. 2011), we 
use the Bouwens et al. (2015) reductions of the HST 
observations over all five CANDELS fields, the version 
7 reduction of the deep CFHT legacy survey observa¬ 


3.1. Catalog Construction and Photometry 

As in previous work (e.g., Bouwens et al. 2007, 2011b, 
2015), we use a modified version of the SExtractor soft¬ 
ware (Bertin & Arnouts 1996) for constructing our HST 
source catalogs that lie at the core of the z ~ 9-10 selec¬ 
tion we perform. SExtractor is run in dual mode, with 
source detection done off the L/i6o-band images and pho¬ 
tometry performed on the V( 306 j Fsi4> J 125 , JHi 4 Ch and 
H 160 images one at a time. Color measurements are 
made in small scalable apertures using Kron-style (1980) 
photometry and a Kron-parameter of 1.2. Fluxes mea¬ 
sured in these small scalable aperture are then corrected 
to total in two steps. In the first step, we multiply each 
of the fluxes by the excess flux seen in the larger scalable 
apertures (Kron parameter of 2.5) for the L/i6o-band over 
that present in smaller scalable apertures. In the second 
step, we correct for the light on the wings of the PSF 
and outside our larger-scalable apertures, based on the 
tabulated encircled energy corrections for point sources 
(Dressel et al. 2012). 

For measurements of the flux in the ground-based ob¬ 
servations or the Spitzer /IRAC observations, we use the 
mophongo software (Labbe et al. 2006, 2010a,b, 2013, 
2015). This software allows to cope with the significant 
amounts of overlap in the light distribution for nearby 
sources. As with other software in the literature with 
similar objectives, MOPHONGO attempts to overcome the 
issue of source confusion by assuming the high-resolution 
HST images (here the Hi6o-band image) provide an ac¬ 
curate model of spatial profile of sources in the ground- 
based/ Spitzer/TR.AG images and that only the normal¬ 
ization of source flux varies from one passband to an¬ 
other. MOPHONGO then varies their individual fluxes to 

11 http://www.cfht.hawaii.edu/Science/CFHTLS 
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Fig. 2.— iTi6o-band images of the CANDELS-UDS, COSMOS, and EGS search data that we used to identify tentative candidate z ~ 9-10 
galaxies. The positions of the candidate z ~ 9-10 sources are indicated by the stars on these mosaics. The numbers adjacent to the stars 
indicate the identity of the tentative z ~ 9-10 candidate (identical the numbering scheme employed in Table [3]). Those stars shaded in grey 
indicate sources that were explicitly preselected for targeted follow-up observations from our 11-orbit z9-CANDELS program, while those 
shaded in yellow were not preselected and only incidentally targeted (lowest two rows in Table Also shown are the regions within these 
fields where deep near-IR observations are available from programs like HUGS (Fontana et al. 2014) or ZFOURGE (Labbe et al. 2015, in 
prep) which cover most but not all of the area over the targeted CANDELS fields. The candidates enclosed in red or purple circles appear 
very likely (>90% confidence) to be at 2 ~ 9 or z ~ 10, respectively, based on the available photometric constraints (obtained with our 
HST follow-up program or archival observations). 
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TABLE 2 

Selection Criteria Utilized in Assembling our z ~ 9 and 2 ~ 10 Samples 


Redshift Selection Criteria 

Sample Preselection for Targeted HST Follow-up After HST Follow-up 


9 


10 

9 


10 


CANDELS-UDS + CANDELS-COSMOS 
(J125 — #160 > 0.5) A (tfieo — [3-6] < 1.4)A (S/N in both +606 and /gi 4 < 2)A 
(rms S/N in Vfeoe and I 8 14 < 1)A (S/N(H 160 )> 5) A (Xjh 140 h 160 ) > 36 ) A 
(^•K,[3.6 ],[4. 5] '' > A pre(^ > 8) > 0.5) A (P post( z > 8) > 0.9) 

idem 


CANDELS-EGS 

(J 125 — #160 > 0.5) A (#160 — [ 3 . 6 ] < 1.4)A (S/N in both +606 and Lsil < 2)A 
(rms S/N in V 606 and 7 8 14 < 1)A (S/N(H 160 )> 5) A (xjh 140 h 160 ) > 6 ) A 
(^K,[3.6],[4.5] ^ A ( Ppost( z > 8) > 0.9) 

idem 


(P(2 > 8) > 0.9)A(8.4 < Zphot < 9.5) 


(P (2 > 8) > 0.9) A (9.5 < z vhot < 11) 
(P(z > 8) > 0.9)A(8.4 < z vhot < 9-5) 


(P(2 > 8) > 0.9) A (9.5 < Zp hot < 11) 


CANDELS-GOODS-North + CANDELS-GOODS-South + ERS 

9 ((Y-dropout criterion from Bouwens+2015) V(Ji 25 — #160 > 0.5)) A (#160 — [3-6] < 1.4)A 

(S/N in both +606 and I B 14 < 2) A (x | 435 T 606 i 775 / 814 , 850 < 4 ) A 
(P(2 > 8) > 0.8)A(8.4 < Zp hot < 9.5) 

10 ((Y-dropout criterion from Bouwens+2015) V(Ji 25 — #160 > 0.5)) A (#160 — [3.6] < 1.4)A 

(S/N in both +606 and / 8 i 4 < 2) A (x| 436 v 606 i 776 I 814 , 860 < 4 )^ 

_ (P(z > 8) > 0.8)A(9.5 < Zphot < 11-0)_ 

Redshift likelihood probability P(z) are computed using our flux meaurements in all photometric bands listed in Table [l] P pr e{z > 8) 
indicates the probability that a source has a redshift greater than 8 before acquiring any follow-up observations, while P p0 st(z > 8) indicates 
the probability that a source has a redshift greater than 8, after obtaining the 1-orbit of follow-up HST observations (assuming the follow-up 
observations yielded a measured flux of 0d=12 njy in the Yios-band filter). 


TABLE 3 

z ~ 9-10 Candidate Galaxies over the CANDELS UDS, COSMOS, and EGS program targeted with our Z9-CANDELS 

FOLLOW-UP PROGRAM. 


ID 

R.A. 

Dec 

#160, AB 

z P hot, P re 

P pre{z > 8) a 

r b 

/ C , P hot, P ost 

P post (z > 8) b 


z 

= 9-10 Candidates Preselected for Targeted Follow-Up Observations with HST 


COS910-0 

10:00:43.16 

02:25:10.5 

26.2+0.1 

9.1 

0.72 

7.8 

0.47 

COS910-I 

10:00:30.34 

02:23:01.6 

26.4+0.2 

9.0 

0.95 

9.0 

0.99 

COS9IO-2 

10:00:14.91 

02:12:10.8 

26.3+0.2 

9.3 

0.74 

9.3 

0.37 

COS910-3 

10:00:27.98 

02:11:49.5 

25.9+0.1 

9.2 

0.63 

2.3 

0.27 

UDS910-0 

02:17:55.50 

-05:11:41.3 

26.4+0.2 

8.8 

0.72 

1.7 

0.09 

UDS910-1 

02:17:21.96 

-05:08:14.7 

26.6+0.2 

8.7 

0.74 

8.6 

0.74 

EGS910-0 

14:20:23.47 

53:01:30.5 

26.2+0.1 

9.1 

0.67 

9.1 

0.92 

EGS910-1 

14:20:21.54 

52:57:58.4 

26.6+0.1 

8.9 

0.19+ 

0.4 

0.02 

EGS910-2 

14:20:44.31 

52:58:54.4 

26.7+0.2 

9.6 

0.69 

9.6 

0.71 

EGS910-3 

14:19:45.28 

52:54:42.5 

26.4+0.2 

8.9 

0.64 

9.0 

0.97 

EGS910-4 

14:19:23.59 

52:49:23.4 

26.2+0.2 

9.2 

0.10+ 

1.0 

0.02 

EGS910-5 

14:19:11.08 

52:46:25.7 

25.8+0.1 

9.2 

0.28 

1.8 

0.11 


z 

~ 9-10 Galaxy Candidates Targeted at No Additional Cost (Not Preselected) 0 


EGS910-6 

14:19:13.84 

52:50:44.7 

26.6+0.2 

9.3 

0.40 

7.0 

0.00 

EGS910-7 

14:20:23.72 

53:01:38.3 

26.0+0.1 

— 

-- 

2.4 

0.18 


a Best-fit z > 4 redshift and integrated z > 8 likelihood for source derived from our HST-\-Spitzer/ IRAC+ground-based photometry 
(Table QJ, before obtaining observations from our z9-CANDELS follow-up program. 

b Best-fit redshift and integrated z > 8 likelihood for source derived from our i/^T+S'pzteer/IRAC-bground-based photometry (Table [TJ, 
after obtaining observations from our z9-CANDELS follow-up program. 

c These sources could be fit within the same WFC3/IR tiles, as our primary targets, and hence required no additional HST time to 
investigate. 

' Over the CANDELS EGS field, we selected sources which, if they showed a null detection in the Y 105 band in a 1-orbit integration, could 
be confirmed with >90% probability to lie at z > 8. While these two sources initially only showed a modest probability for being at z > 8, 
their SEDs were nevertheless consistent with lying at z > 8 (particularly if a null detection at 1.05/mi could be confirmed). 
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obtain a good fit. Measurements of the flux for individ¬ 
ual sources is then performed in fixed circular apertures, 
after subtracting the model light profile from neighboring 
sources. 1 . 2 "-diameter, 1 . 8 "-diameter, and 2 "-diameter 
apertures are used for the ground-based photometry, 
Spitzer/IRAG, and Spitzer /IRAC photometry over all 
of our fields, the CANDELS GOODS-North+GOODS- 
South fields, and the CANDELS UDS/COSMOS/EGS 
fields. The measured fluxes are then corrected to to¬ 
tal based on the model profile for individual sources. 
Narrower apertures are used for our Spitzer/IRAG pho¬ 
tometry over the GOODS-North+GOODS-South fields 
to leverage the narrower FWHM of the Spitzer/IRAG 
PSF in the Labbe et al. (2015) reductions. 

3.2. Selection of Bright z ~ 9 -10 Candidates over the 
CANDELS-UDS and CANDELS-COSMOS Fields 
3.2.1. Selection Criteria 

In searching for candidate galaxies at 2 ~ 9-10, we 
suppose that these galaxies have very similar colors and 
SEDs as galaxies at slightly lower redshifts. Specifically, 
we would expect these sources to show a sharp spectral 
break at 1216A, due to strong absorption from the neu¬ 
tral hydrogen forest, and to exhibit a blue I/U-continuum 
redward of the break. 

For star-forming galaxies at redshifts z ~ 8.4 and 
higher, the Lyman-break will already have redshifted a 
significant way through the Ji 25 -band, yielding moder¬ 
ately red J 125 — ifi 6 o colors. As a result, the selection of 
all sources with red J 125 — -H 16 O colors should allow us to 
identify the bulk of star-forming galaxies from z ~ 8.7 
to z ~ 11 (particularly if those galaxies are not substan¬ 
tially dust obscured). 

Here we search for candidate z > 8.4 galaxies over 
the CANDELS-UDS and COSMOS fields using a J 125 — 
H 160 > 0.5 criterion. Star-forming galaxies with a UV- 
continuum slope /3 of — 1.6 (typical of luminous galaxies 
at z = 4-7) would have J 125 — i /160 color of ~ 0.5 at 
z = 8.7, but the lower-redshift limit for our selection will 
depend on the intrinsic colors of individual galaxies and 
also can be affected by observational noise. 

In addition to our J 125 — Hi 60 criterion, we also re¬ 
quire that sources be undetected (< 2 cr) in the Veso 6 or 
/ 814 -bands. Sources where the root mean square S/N in 
the 1/606 and Is 14 bands is greater than 1 are excluded. 
In addition to these non-detection requirements on the 
HST optical data, we also require that sources remain un¬ 
detected (< 2.5 ct) in an inverse-variance-weighted mean 
stack of the ground-based optical data. 

We also demand that sources show Z/i6o — [3.6] colors 
bluer than 1.4 mag to exclude intrinsically-red or old 
z ~ 2 galaxies from our samples, similar to the criteria 
that Oesch et al. (2014) or Bouwens et al. (2015) apply. 
This particular color cut corresponds to a I/U-continuum 
slope j3 of 0.0 (where f\ oc A^), which is approximately 
as red as bright galaxies area observed to be at z ~ 6-8 
(e.g., Bouwens et al. 2012, 2014a; Wilkins et al. 2011; 
Finkelstein et al. 2012; Rogers et al. 2014). 

We require that all selected z ~ 9-10 candidates show 
strong evidence of corresponding to real sources. We 
therefore require that (1) sources be detected in the Z /160 
band at >5cr significance, (2) the root mean square detec¬ 
tion significance of sources in the J //140 and Z/i6o-band 


images be at least 6, and (3) the root mean square de¬ 
tection significance of sources in the J//i 4 o, [3.6], [4.5], 
and iv-bands be at least 2er. 

Finally, in the last step, we compute the redsliift like¬ 
lihood distribution for each candidate source using the 
EAZY photometric redshift code (Brammer et al. 2008) 
based on the photometry we have available for sources, 
the standard EAZY_vl.O template set, and a flat prior. 
We supplemented the standard EAZY.vl.O template set 
with SED templates from the Galaxy Evolutionary Syn¬ 
thesis Models (GALEV: Kotulla et al. 2009). Nebular 
continuum and emission lines were added to the later 
templates using the Anders & Fritze-v. Alvensleben 
(2003) prescription, a 0 -2Zq metallicity, and a rest-frame 
EW for Ha of 1300A (which appears to be appropriate 
for z ~ 6-7 (Smit et al. 2014, 2015; Roberts-Borsani et 
al. 2015). 

The photometry utilized for constraining the like¬ 
lihood distributions for individual sources included 
the HST V 606 / 8 i 4 ^i 25 -///i 4 o//i 60 + Subaru-SuprimeCam 
BgVriz + CFHT/Megacam ugriyz + UltraVISTA 
YJHK S + ZFOURGE A J 2 J 3 //,//, + Spitzer/IRAG 
3.6/j,m+4.5^m S-CANDELS data sets for the CAN- 
DELS COSMOS held, HST VeoehuJ^JH^oH^o + 
Subaru-SuprimeCam BVriz + CFHT/Megacam 11 + 
UKIRT/WFCAM JHK S + ZFOURGE J 1 J. 2 J :i H s H l 
+ VLT/HAWKI/HUGS YI\ S data sets for the CAN- 
DELS UDS held, and the HST V 606 1 S u J 125 JH U 0 H 160 
+ CFHT/Megacam ugriyz + CFHT/WIRCam K s + 
Spitzer/ IRAC 3.6^m+4.5/mi data sets for the CAN- 
DELS EGS held. The depths of these observations are 
provided in Tableland their areal coverage is illustrated 
in Figure [2] 

Sources that satished our aforementioned criteria, 
which showed a >50% probability of being at z > 8, and 
which could be confirmed to be a >90% likelihood can¬ 
didate with a single orbit of follow-up observations (sup¬ 
posing sources are measured to have a hux of Oil2 nJy 
in the Yios band) made it into our final preselection of 
candidate z ~ 9-10 galaxies (to be targeted with follow¬ 
up observations). In computing the posterior probability 
that a source has a redshift z > 8 or z < 8, we adopt a 
hat prior on the redshift. 

Table [2] provides a convenient compilation of all the 
selection criteria we employed in preselecting candidate 
z ~ 9-10 galaxies to follow-up with targeted observa¬ 
tions. 

3.2.2. UDS+ COSMOS Results 

Applying the selection criteria from the previous sec¬ 
tion to our source catalogs over the CANDELS-UDS and 
CANDELS-COSMOS fields, we found five sources which 
satisfied all of the criteria. A list of all 6 sources satis¬ 
fying these criteria are included in Table [3] and similar 
candidates from the CANDELS-EGS held. 

The observed spectral energy distributions for these 
hve candidate z ~ 9-10 galaxies are presented in Ap¬ 
pendix A (Figure [H]), along with SED hts to a model 
z > 8 galaxy and a model z < 3 galaxy. Also shown 
on this figure is the redshift likelihood distribution ( solid 
black line) based on the photometry we have available for 
each candidate in the ~20 different wavelength channels 
(HST + Spitzer/IRAG + ground-based observations). In 
addition, this figure presents the redshift likelihood dis- 
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Fig. 3.— HST + Spitzer /IRAC images for 5 candidate 2 ~ 9-10 galaxies which were confirmed as probable 2 : > 9 galaxies (or partially 
confirmed in the case of EGS910-2) using HST follow-up observations with our z9-CANDELS program. Fits to the SEDs of these sources 
and the estimated redshift likelihood distributions are presented in Figure EO 


TABLE 4 

Photometrically-Confirmed 2 ~ 9-10 Galaxies over the CANDELS Fields 


ID 

R.A. 

Dec 

Hl60,AB 

z phot° 

P (2 > 8 ) 

Ref a 



2 

~ 9 Sample 




Two-Part Search Strategy (Preselection + hollow-up: §3, §4): 


9 O " 1 ” 0 ' 4 

y ' u -0.5 



COS910-1 

10:00:30.34 

02:23:01.6 

26.4±0.2 

0.99 


EGS910-0 

14:20:23.47 

53:01:30.5 

26.2±0.1 

Q 1+ 0 - 3 

y,± -0.4 

0.92 


EGS910-3 

14:19:45.28 

52:54:42.5 

26.4±0.2 

9 Q + 0 - 5 
y ' u -0.7 

0.97 


UDS910-1 C 

02:17:21.96 

-05:08:14.7 

26.6±0.2 

0 c+ 0.6 

®- b -0.5 

0.74 


Direct Search 

Strategy for 2 

> 8.4 Galaxies (§5): 





GS-z9-l 

03:32:32.05 

-27:50:41.7 

26.6±0.2 

9.3±0.5 

0.9992 

[1]. PI 

GS-z9-2 

03:32:37.79 

-27:42:34.4 

26.9 

8-9to:s 

0.83 

[2] 

GS-z9-3 

03:32:34.99 

-27:49:21.6 

26.9 

0 O+0.3 

0.3 

0.95 

[2], [3] 

GS-z9-4 

03:33:07.58 

-27:50:55.0 

26.8 


0.97 

[2], [3] 

GS-z9-5 

03:32:39.96 

-27:42:01.9 

26.4 

8.7i°;8 

0.55 

[2] 

GN-z9-l 

12:36:52.25 

62:18:42.4 

26.6±0.1 

9.2±0.3 

>0.9999 

[1], [2] 



2 

~ 10 Sample 




Two-Part Search Strategy (Preselection + Follow-up: §3, §4): 


9 0+°- 5 

0.5 



EGS910-2 C 

14:20:44.31 

52:58:54.4 

26.7±0.2 

0.71 


Direct Search Strategy for 2 

> 8.4 Galaxies (§5): 





GN-zl0-l d 

12:36:25.46 

62:14:31.4 

26.0±0.1 

ll.liO.l 

>0.9999 

[1], [2], [4], [5] 

GN-zlO-2 

12:37:22.74 

62:14:22.4 

26.8±0.1 

9.9±0.3 

0.9994 

[1], [2] 

GN-zlO-3 

12:36:04.09 

62:14:29.6 

26.8±0.2 

9.5±0.4 

0.9981 

[1], [2] 

GS-zlO-1 

03:32:26.97 

-27:46:28.3 

26.9±0.2 

9.9±0.5 

0.9988 

[1 , PI 


a References: [1] Oesch et al. 2014, [2] Bouwens et al. 2015, [3] McLure et al. 2013, [4] Oesch et al. 2016, [5] Bouwens et al. 2010 
b lcr uncertainties are computed based on the z > 4 likelihood distributions. 

c This candidate could only be partially confirmed, given the limited orbit allocation to our HST program. 

d This source is now spectroscopically confirmed to lie at 2 = 11.1 (Oesch et al. 2016), but broadly lies within our 2 ~ 10 selection window. 
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Fig. 4. — {left) Best-fit SED models to the observed i/5T+5pi£zer/IRAC-|-ground-based photometry of five candidate z ~ 9-10 galaxies 
(COS910-0, UDS910-1, EGS910-0, EGS910-2, EGS910-3) that have been photometrically confirmed (or partially confirmed in the case 
of UDS910-1 and EGS910-2) by observations from the z9-CANDELS follow-up program. Red solid circles, la error bars, and la limits 
are from the HST or Spitzer /IRAC observations, while the black solid circles, la error bars, and la limits are from the ground-based 
observations. The solid blue line shows the best-fitting SED for a z > 6 galaxy, while the grey line shows the best-fitting SED for a 
z < 6 galaxy, {right) Redshift likelihood distribution for these z ~ 9-10 candidates incorporating both our follow-up observations and the 
HST-\- Spitzer / IRAC+ground-based observations that were used in the pre-selection {solid lines). 
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tribution we would expect, assuming these candidates are 
not detected in the single orbit of follow-up Eios-band 
observations from the z9-CANDELS program. 

Postage stamp images of these six candidates are also 
presented in Appendix A ('Figure HOl). As should be ob¬ 
vious from this figure, all six of the present z ~ 9-10 
candidates show clear detections in the Hi 6 o-band, as 
well as significant ~ 2-3 a detections in Ji 25 -band and 
Ji?i 40 -band observations (where available), as well as in 
the S-CANDELS Spitzer/IRAC data. 

All six of these candidates bear a remarkable similarity 
to the first samples of particularly luminous z ~ 9-10 
galaxies identified by Oescli et al. (2014) in terms of their 
very blue Hi eo — [3.6] colors (see also Wilkins et al. 2016), 
red [3.6] — [4.5] colors, and observed sizes (Holwerda et 
al. 2015). 

3.3. Selection of z ~ 9 -10 Candidates over the 
CANDELS-EGS 

3.3.1. Selection Criteria 

The selection of candidate z ~ 9-10 galaxies over the 
CANDELS-EGS field is even more challenging than over 
the CANDELS-UDS and CANDELS-COSMOS fields 
due to the lack of deep observations at 1.05/zm over the 
CANDELS-EGS field. E-band observations (at 1.05/zm) 
play a crucial role in excluding the possibility that 
sources can correspond to slightly reddened star-forming 
galaxies at z ~ 7.5-8.5 or correspond to passive or red¬ 
dened galaxies at much lower redshifts. 

In selecting z ~ 9-10 candidates over the CANDELS- 
EGS field, we therefore adopted almost identical crite¬ 
ria as over the CANDELS-COSMOS or CANDELS-UDS 
fields, with one exception. Instead of requiring sources 
to have a >50% probability of corresponding to a 2 > 8 
galaxy, we required that sources be capable of confirma¬ 
tion with a single orbit of HST observations at 1.05/zm. 
For the purposes of selection, we take confirmation to 
correspond to the source having >90% likelihood of be¬ 
ing at z > 8 after adding a flux constraint of 0±12 nJy 
to the observed SED at 1.05/xm (though we obtained 
follow-up observations in JH\^q for the one case where 
the J 125 — -ffi 60 color was > 1 . 2 ). 

3.3.2. EGS Results 

Applying the selection criteria from the previous sec¬ 
tion to our source catalogs over the CANDELS-EGS 
field, we found six additional sources which satisfied all 
of the criteria (Table [3]). 

The observed spectral energy distributions for these six 
candidate z ~ 9-10 galaxies are presented in Appendix A 
('Figure fill), along with SED fits to a model z> & galaxy 
and a model z < 6 galaxy. Postage stamp images of these 
six candidate z ~ 9-10 galaxies are also provided. 

The most promising z ~ 9-10 candidates we identi¬ 
fied over the CANDELS-EGS field were the EGS910-0, 
EGS910-2, and EGS910-3. All 3 sources show evidence 
for a sharp break break at 1 . 2 /xm, as well as a blue spec¬ 
tral slope redward of the break. The other candidates 
also show evidence for a strong spectral break at 1 . 2 /xm 
and a blue spectral slope redward of the break, but also 
show possible flux in ~l -2 passbands blueward of the 
break. Until observations from our z9-CANDELS follow¬ 
up program became available on these candidates, it was 


not possible to determine whether they are more likely 
to correspond to bona-fide z ~ 9-10 galaxies or z ~ 1-3 
interlopers. 

4. NATURE OF THE TARGETED Z ~ 9-10 CANDIDATES 

HST observations are now available over all 12 can¬ 
didate z ~ 9-10 galaxies targeted by our z9-CANDELS 
program. These observations allow us to make a fairly 
definitive assessment of the nature of these candidate 
z ~ 9-10 based on the flux we measure for these can¬ 
didates at 1.05/im. One orbit of E 05 -band observa¬ 
tions have already been obtained for eight candidates tar¬ 
geted by our program COS910-0, COS910-1, COS910-3, 
UDS910-0, UDS910-1, EGS910-0, EGS910-1, EGS910- 
3, and EGS910-4. Slightly shallower observations (i.e., 
1/3 and 2/3 of an orbit) in the Eios-band were acquired 
on the candidates EGS910-5 and COS910-2, due to the 
greater brightness of the former candidate and the util¬ 
ity of an additional 1/3 orbit Ji/ 140 -band observations to 
investigate the nature of the potential z ~ 10 candidate 
galaxies EGS910-2 and COS910-3. 

Eo 5 _ band images for these candidates are presented in 
either Figures 151 or 1121 from Appendix B, in conjunction 
with images of these candidates at other wavelengths. 
Figures 0] and [13] from Appendix B show the observed 
SEDs for the targetted z ~ 9-10 candidates in CANDELS 
program. 

The present observations confirm photometrically 5 of 
the first 12 z ~ 9-10 candidates targeted by our pro¬ 
gram. Two of these five confirmations are only partial 
confirmations (EGS910-2 and UDS910-1: as more obser¬ 
vations are needed for these candidates to >90% secure). 
Detailed remarks on the confirmed z ~ 9-10 candidates 
can be found here: 

COS910-1: COS910-1 is not detected (< lcr) in the 
Elos-band follow-up observations at 1.05/zm. A detailed 
fit to its SED suggests that it is actually a star-forming 
galaxy at z = 9.1, with <0.7% probability of it corre¬ 
sponding to a 2 < 8 galaxy. 

UDS910-1: UDS910-1 is not detected at (<lcr) in the 
E. 05 -band follow-up observations we obtained at 1.05/xm. 
Rederiving the redshift likelihood distribution using the 
new flux information in the Eios-band, we compute a 
best-fit photometric redshift of 8 . 6 , with a 4% and 24% 
probability of corresponding to a z < 7 and z < 8 source, 
respectively. 

EGS910-0: EGS910-0 is not detected (<1ct) in the Eios- 
band follow-up observations at 1.05/zm. Rederiving the 
redshift likelihood distribution using the new flux infor¬ 
mation in the Eios-band, we compute a best-fit photo¬ 
metric redshift of 9.1, with only a 4% probability of cor¬ 
responding to a z < 8 source. 

EGS910-2: Follow-up of EGS910-2 in the J-HEo-band 
shows a clear 2.6er detection of the source and which is 
in excellent agreement with the expected flux given a 
model redshift of z ~ 9.6 for the source. Nevertheless, 
the source is sufficiently faint that the redshift likelihood 
distribution shows a 29% likelihood of the source being 
at 2 < 8 . Deeper follow-up observations at 1.05/um will 
be required to rule out the z < 8 solution. 

EGS910-3: EGS910-3 shows no detection (<lcr) in the 
Eios-band follow-up observations we obtained at 1.05/rm. 
Rederiving the redshift likelihood distribution using the 
new flux information in the Eios-band, we compute a 
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Fig. 5.— (left) Best-fit SED models to the observed FT5TH-Spi£zer/IRAC+ground-based photometry of three candidate z ~ 9 galaxies 
(GS-z9-2, GS-z9-3, GS-z9-4) that satisfied our criteria for selection. This figure also includes another candidate z ~ 9 galaxy (GS-z9-5), 
but appears to nevertheless be a high-redshift galaxy. Spitzer /IRAC observations are not shown for GS-z9-5 since it is nearby a very 
bright star and under a diffraction spike. These sources were identified in a separate search over the extended GOODS-South area (ERS, 
CANDELS GOODS-South, HUDF09-1, HUDF09-2: see §5.1). The points and lines are otherwise as in Figure[4] (right) Redshift likelihood 
distribution for these z ~ 9 candidates (solid lines). 
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Fig. 6. — HST + Spitzer /IRAC images (6 ;/ x 6 /; ) of three candidate z ~ 9 galaxies (GS-z9-2, GS-z9-3, GS-z9-4) that satisfied our criteria 
for selection and another candidate z ~ 9 galaxy (GS-z9-5) that did not satisfy these criteria (but appears nevertheless to be a high-redshift 
source). Spitzer /IRAC observations are not shown for GS-z9-5 since it is nearby a very bright star and under a diffraction spike. We 
have identified these candidates in a separate search over the extended GOODS-South area (ERS, CANDELS GOODS-South, HUDF09-1, 
HUDF09-2: see §5.1). 
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best-fit photometric redsliift of 9.0, with a 3% probabil¬ 
ity of corresponding to a z < 8 source. 

The 3 j> 8.4 candidates confirmed - and 2 candidates 
partially confirmed - by our z9-CANDELS program are 
compiled for convenience in Table 0] We will also include 
in this table some additional candidates we identify in §5 
(and also as identified by Oesch et al. 2014 and Bouwens 
et al. 2015). 

Our overall confirmation rate is 42% (5/12) for sources 
preselected by our criteria. We achieve an even higher 
56% success rate targeting those sources from our selec¬ 
tion which are high-probability (>50%) 2 > 8 galaxies, 
before our follow-up observations. While imperfect, this 
program is very efficient, supplementing some 270 orbits 
of HST time and hundreds of hours of Spitzer time with 
only 11 orbits of additional HST time. By contrast, the 
CANDELS + ERS programs over the GOODS-North + 
South cost some ~500 orbits, and we identified only 9 
candidates in those data, or 0.02 z ~ 9-10 candidate per 
invested orbit. 

Detailed remarks on the 2 ~ 9-10 candidate galaxies 
which were not confirmed by our follow-up program are 
included in Appendix A. 

5. COMPLETING THE CENSUS OF CANDIDATE Z ~ 9 
GALAXIES OVER THE CANDELS GOODS-NORTH, 
GOODS-SOUTH, AND ERS FIELDS 

We can obtain the best constraints on the volume den¬ 
sity of bright z ~ 9-10 galaxy candidates by not sim¬ 
ply considering a search over the CANDELS UDS, COS¬ 
MOS, and EGS fields as we did in the previous sections, 
but also considering a search for similar sources over the 
GOODS-North and GOODS-South fields. 

5.1. Criteria for Identifying z ~ 8.5 -9.0 Galaxies 

The purpose of the present section is to obtain a com¬ 
plete census of the bright z ~ 8.5-11 galaxy candidates 
over the CANDELS GOODS-North+GOODS-South + 
ERS fields. 

In Oesch et al. (2014) and Bouwens et al. (2015), 
we already conducted a significant search for galaxies 
in this redshift range, by looking for sources with red 
<7i25 — Hiqo > 0.5 colors and blue H i 60 — [3.6] < 1.4 
colors. However, such a selection is only sensitive to 
galaxies with redshifts z > 9 and can suffer significant 
incompleteness at z < 9. 

Here we extend the search from Oesch et al. (2014) and 
Bouwens et al. (2015) to also consider sources with red- 
shifts z > 8.4. We select these sources by considering all 
those sources which satisfy the z ~ 8 color-color criteria 
of Bouwens et al. (2015), deriving photometric redshifts 
for all such sources using the EAZY photometric redshift 
code (Brammer et al. 2008), and including those sources 
where the most likely redshift is greater than 8.4. 

The photometry we consider in deriv¬ 
ing the redshift likelihood contours are the 
Bouwens et al. (2015) reductions of the HST 
I?435^606*775^814~850l / 098^ / 105 J 125 J -^140^160 data, 

the Labbe et al. (2015) reductions of essentially all 
Spitzer/TRAC observations over the GOODS-North and 
South fields, and the Bouwens et al. (2015) reductions 
of the HUGS HAWK-I A' s -band observations. 

Briefly, the Bouwens et al. (2015) selection criteria for 


identifying z ~ 8 sources is 

(Ul 05 — J\ 2 $ > 0.45) A (J 125 — #160 < 0.5) A 
(Ulos ~ Ji25 > 0.75( J 125 — #L6o) + 0.525) 

for sources over the CANDELS GOODS-North + 
GOODS-South fields and 

(^098 ~ Jl 25 > 1-3) A (Ji25 — #160 < 0.5) A 
(^098 — J125 > 0.75( J 125 — #160 ) + 1.3) 

for sources over the ~40 arcmin 2 ERS field. Sources are 
required to be detected at 6cr in a y 2 stack of the H 160 - 
band or JH^q + #160 band observations redward of the 
break (in a fixed O^O^-diameter aperture). 

To ensure contamination is kept to a minimum, an 
optical “x 2 ” is computed for each candidate source 
(Bouwens et al. 2011b) based on the flux in the 
1 ^ 435 ^ 606 ^ 775 ^ 814 -^ 850 -dand observations. xl P t ^ taken to 
equal EjSGN(/ i )(/ i /cr i ) 2 where fi is the flux in band i in 
a consistent aperture, cq is the uncertainty in this flux, 
and SGN(/j) is equal to 1 if fi > 0 and —1 if fj < 0. Any 
candidate with a measured \ op t in excess of 4 is excluded 
from our selections. 

We only search for z > 8.4 sources over the GOODS- 
North and GOODS-South fields brightward of #i 6 o,a,b = 
27 mag to ensure that we have strong constraints on 
the nature of the selected sources to the limit of our 
search. The effective depth of our 2 ~ 9-10 search over 
the CANDELS-UDS, COSMOS, and EGS fields is also 
approximately ~27 mag, so the effective depth of our 
search is similar across all five CANDELS fields that we 
utilize. 

5.2. Selection Results 

Using the selection criteria from the previous section, 
we identify three high-probability (>90% confidence) and 
one moderate probability (~50% confidence) z ~ 8 .4-9.0 
galaxies over the ERS, CANDELS GOODS-South, and 
CANDELS GOODS-South fields. 

The 17i6o-band magnitudes of the z ~ 8.4-9.0 galaxies 
we have selected range from 26.4 and 26.9, similar to that 
found for our z ~ 9-10 sample over the CANDELS-UDS, 
COSMOS, and EGS fields. We have included the four 
new z ~ 8 .4-9.0 candidates in Table 01 along with other 
high-probability z ~ 8.4-11 identified here. Fits to the 
observed SEDs for our new z ~ 8 .4-9.0 candidates over 
these fields are shown in Figured Postage stamp images 
of the candidates are provided in Figure EE 

5.3. Criteria for Identifying z > 9.0 Galaxies 

As performed by Oesch et al. (2014) and Bouwens et 
al. (2015), we also include sources with J i2 5 — H i 60 > 
0.5, H 160 — [3.6] < 1.4 colors. Our selection criteria for 
identifying these sources are essentially identical to that 
utilized by Bouwens et al. (2015: see also Oesch et al. 
2014) except that a J 125 — #160 > 0.5 color criterion is 
utilized. 

We identify exactly the same set of sources as Oesch 
et al. (2014) identify using the above criteria. A com¬ 
pilation of these sources and other z ~ 8 .4-9.0 sources 
identified over the CANDELS GOODS-North, GOODS- 
South, and ERS fields is provided in Table [4] 
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Fig . 7 . — Simple binned determinations of the UV LF for lumi¬ 
nous galaxies at z = 9 (upper panel) and z = 10 (lower panel), lcr 
upper limits on the volume density of z ~ 9 and z ~ 10 galaxies 
are included at 22 mag. The shaded hatched red region in¬ 
dicates the volume densities (at a given Mjjy) preferred at 68% 
confidence by this analysis (see Table [6]). To put these constraints 
on the bright end of the UV LF in context, we also include determi¬ 
nations of the z ~ 9 and z ~ 10 UV LFs at lower luminosities from 
Zheng et al. (2012: solid blue square), Oesch et al. (2013: open red 
square), McLure et al. (2013: open blue triangles), and McLeod 
et al. (2015: open blue squares). The lightly-shaded red region 
shows the constraints on the z ~ 9 LF as derived by Bouwens 
et al. (2014b) using a search for z ~ 9 galaxies over the CLASH 
program (Postman et al. 2012). The overplotted line shows an ex¬ 
trapolation of the Bouwens et al. (2015) LF results to £ ~ 9 and 
z rv* 10.2 based on the fitting formula provided in §7.1. 


6. IMPACT OF GRAVITATIONAL LENSING FROM 
FOREGROUND GALAXIES ON THESE RESULTS 

From previous work (Wyithe et al. 2011; Barone- 
Nugent et al. 2015; Mason et al. 2015; Fialkov & Loeb 
2015), it is well known that gravitational lensing from 
foreground galaxies can have a particularly significant 
effect in enhancing the surface density of bright z > 6 
galaxies on the sky. This is especially true for the bright¬ 
est sources due to the intrinsic rarity and the large path 
length available for lensing by foreground sources. 

Given this phenomenon, it has become increasingly 
common for researchers searching for the brightest z ~ 6- 
10 galaxies to look for possible evidence of lensing am¬ 
plification (Oesch et al. 2014; Bowler et al. 2014, 2015; 
Zitrin et al. 2015; Roberts-Borsani et al. 2016). While 


TABLE 5 

Binned Determination of the rest-frame UV LF at z ~ 9 
AND 2 ~ 10 . 


Ati 600 ,AB a <t>k (10 3 Mpc 3 mag 1 ) 


z ~ 9 galaxies 


- 21.94 


< 0 . 0024 b 

- 21.14 


0 0044" 1 ” 0 ' 0042 
u,uu ^ -0.0024 

- 20.34 

z ~ 10 galaxies 

0 0322~ l ~°‘° 217 

u . uo ^_ 0 . 0138 

- 22.05 


< 0 . 0017 b 

- 21.25 


0 . 0009 ± 8 ; 8 ™ 

- 20.45 


0 . 0180 i °; 3 M 8 


a Derived at a rest-frame wavelength of 1600A. 
b lcr upper limit. 

there are a number of cases where such magnification 
boosts may be present (e.g., Barone-Nugent et al. 2015; 
Roberts-Borsani et al. 2016), the fraction of lensed 
sources among bright samples still does not appear to 
be particularly high (Bowler et al. 2015). 

As in the above work, we explicitly check our compi¬ 
lation of bright z ~ 9-10 galaxy candidates from these 
fields for evidence of gravitational lensing. For conve¬ 
nience, we use the Skelton et al. (2014) catalogs provid¬ 
ing radii and stellar mass estimates for all sources over 
the CANDELS areas we have searched. The Skelton et 
al. (2014) catalogs use the diverse multi-wavelength data 
over the CANDELS fields, including HST optical, near- 
infrared, Spitzer/IRAC, and ground-based observations, 
to provide flux measurements of a wide wavelength range 
and then use these flux measurements to estimate the 
redsliifts and stellar masses. 

As in Roberts-Borsani et al. (2016), we model galaxies 
in our bright z ~ 9-10 sample as singular isothermal 
spheres, and we use the measured half-light radius and 
inferred stellar mass to derive a velocity dispersion 
estimates for individual galaxies in these samples. We 
found only two examples of galaxies whose measured 
fluxes appear likely to be slightly boosted (>0.1 mag) 
by lensing amplification: 

EGS910-3: There is a foreground galaxy at z ~ 1.9 with 
estimated stellar mass of a 10 1o ' 32 -Mq that lies within 
1.9 arcsec of this source. Based on the velocity dispersion 
we estimate for this source of ~220 km/s, we compute a 
magnification boost of 0.25 mag for this source. 
GN-zlO-2: This source is estimated to be boosted by 
0.11 mag by a 10 1064 Mq galaxy with a spectroscopic 
redsliift of z = 1.02 (Barger et al. 2008) that lies within 
4.0" of the targeted source. This source was previously 
flagged as being slightly lensed by Oesch et al. (2014). 

7. IMPLICATIONS OF OUR SEARCH RESULTS 

7.1. Constraints on the UV Luminosity Functions at 
z = 9 and z = 10 

In this section, we utilize the combined sample of z ~ 
9-10 candidates over the CANDELS-UDS, CANDELS- 
COSMOS, and CANDELS-EGS fields and similar z ~ 
9-10 candidates over CANDELS GOODS-North and 
GOODS-South fields (§5) to quantify the UV LFs at 
0 ~ 9 and 2 ~ 10. Table 0] provides a compilation of the 
relevant sources for our determination of the LF. 

As in our recent paper on the z ~ 4-10 LFs, we use 
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Fig. 8 .— Present determinations of the bright end of the UV 
LF at 2 : ~ 9 and z ~ 10 using all five CANDELS fields (orange 
and dark purple circles, respectively, with lcr uncertainties), lcr 
upper limits are as in Figure [7] For context, we also include fainter 
determinations of the UV LFs at z ~ 9 and z ~ 10 from Oesch et al. 
(2013) and Bouwens et al. (2015: small circles and solid line). The 
dashed lines indicate extrapolations of the Bouwens et al. (2015) 
LF relations to z ~ 9 and z ~ 10.2 (see §7.1). The 2 : ~ 4, 2 : ~ 5, 
z ~ 6, z ~ 7, and 2 : ~ 8 LF determinations from Bouwens et al. 
(2015) are also shown. 


the results from these simulations to derive the selection 
volumes needed to relate the UV LF function cj)(M ) to 
the observed surface density of sources on the sky. For¬ 
mally, we write the UV LF in stepwise format <pj as 
T,(j>jW(M — Mj ) where j is an index running over the 
magnitude bins, where Mj corresponds to the absolute 
magnitude at the center of each bin, where 

0, x < —0.4 

W{x) = 1, -0.4 < x < 0.4 (1) 

0, x > 0.4. 


and where x gives the position within a magnitude bin. 
We take the width of magnitude bins to be 0.8 mag (e.g. 
versus the 0.5-mag used by Bouwens et al. 2015), given 
the limited number of bright z ~ 9 and z ~ 10 galaxies. 

We then look for the derived LF (pj that yields the 
observed surface density of z ~ 9 and 2 ~ 10 galaxies on 
the sky with maximum probability C. As in Bouwens et 
al. (2015), the likelihood C is computed as 

C = TlfieidUiPimi) ( 2 ) 

where the above products runs over the different search 
fields and magnitude interval i used in the LF determi¬ 
nations, and where p(rrii) is probability of identifying a 
certain number of sources in magnitude interval i in a 
given search field. 

For simplicity (and given the small numbers in each of 
our samples: see the discussion in §4 of Bouwens et al. 
2008), we ignore field-to-field variance in deriving the LF 
results and compute the likelihood that our survey fields 
show a certain number of sources assuming Poissonian 
statistics. We therefore compute p{mi) as follows: 


p{rrii ) 


= e 


-N e: 


(A^exp.i) 


Nobs 


(A obs,i,j ) 


(3) 


where N 0 b St i is the observed number of sources in search 


TABLE 6 

68% CONFIDENCE REGIONS ON THE VOLUME DENSITY OF GALAXIES 
AT 2 ~ 9 AND z ~ 10 VS. M uv 


-^1600 ,AB a 

Volume Density (10 
Lower Bound b 

" 3 Mpc 3 mag x ) 

Upper Bound b 

-21.84 

2 : rsj 9 galaxies 
0.0006 

0.0024 

-21.64 

0.0010 

0.0032 

-21.44 

0.0016 

0.0044 

-21.24 

0.0024 

0.0060 

-21.04 

0.0039 

0.0084 

-20.84 

0.0059 

0.0120 

-20.64 

0.0087 

0.0174 

-20.44 

0.0126 

0.0262 

-20.24 

0.0176 

0.0402 

-20.04 

0.0250 

0.0621 

-21.95 

2 : ~ 10 galaxies 
0.0004 

0.0018 

-21.75 

0.0006 

0.0022 

-21.55 

0.0009 

0.0028 

-21.35 

0.0012 

0.0036 

-21.15 

0.0017 

0.0048 

-20.95 

0.0023 

0.0066 

-20.75 

0.0030 

0.0093 

-20.55 

0.0039 

0.0135 

-20.35 

0.0049 

0.0200 

-20.15 

0.0060 

0.0299 


a Derived at a rest-frame wavelength of 1600A. 
b 68 % Confidence Region 


field and magnitude interval i, where N exp j is the ex¬ 
pected number of sources in a search field and magnitude 
interval i. The expected number of sources in a search 
field N expectedii is computed as 

A expected, i — LI j ( 4 ) 

where is the effective volume over which one could 
expect to find a source of absolute magnitude j in the 
observed magnitude interval i. 

The selection volumes Vij are estimated using almost 
an identical procedure to that in Bouwens et al. (2015). 
Specifically, we constructed catalogs with mock sources 
spanning the entire redsliift range 2 ~ 7.5 to 2 ~ 12. 
To ensure that sources had reasonable sizes and mor¬ 
phologies, we randomly selected similar luminosity z ~ 4 
galaxies from the Hubble Ultra Deep Field (Beckwith et 
al. 2006; Illingworth et al. 2013) to use as a template to 
modeling the two-dimensional pixel-by-pixel profiles of 
individual sources. The size of the model sources were 
assumed to scale with redshift as (1 + z)~ 12 to match 
the size scaling observed for sources with fixed luminos¬ 
ity from 2 ~ 10 to 2 ~ 2 (Bouwens et al. 2015; Oesch et 
al. 2010; Ono et al. 2013; Holwerda et al. 2014; Kawa- 
rnata et al. 2015; Shibuya et al. 2015). UV continuum 
slopes of sources were assumed to have a mean value of 
— 1.8, consistent with that measured at high luminosities 
at 2 ~ 5-8 (Bouwens et al. 2012, 2014; Finkelstein et 
al. 2012; Willott et al. 2013; Rogers et al. 2014), with a 
dispersion of 0.3 (Bouwens et al. 2012; Castellano et al. 
2012 ). 

We generate simulated images of each source in all 
HST , ground-based, and Spitzer /IRAC wavelength chan¬ 
nels. Artificial images of individual sources in the 
ground-based and Spitzer /IRAC channels are produced 
by convolving the simulated HST images with the PSF- 
matching kernels we derive with MOPHONGO (Labbe et 
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al. 2013). These images are then added to sections of 
the CANDELS UDS, COSMOS, EGS, GOODS-North, 
and GOODS-South, and ERS fields, catalogs are con¬ 
structed, and sources are selected using exactly the same 
procedures as we apply to the real observations. We in¬ 
clude both the criteria used for our pre-selection and our 
confirmation criteria (i.e., P(z > 8) > 0.9) in computing 
the selection volume. We implement these criteria in an 
identical way as they are applied to the observations. 

For example, to be included in our selection volume 
estimates, simulated sources are determined to be pre¬ 
selected using the criteria we describe in §3.2.1 or §3.3.1. 
For simulated sources within the CANDELS-UDS and 
CANDELS-COSMOS data sets, this means that their 
cumulative probability of lying at z > 8 must be greater 
than 50%, before the addition of any Yios-band data. 
In addition, simulated sources (over the CANDELS- 
UDS/COSMOS/EGS fields) must show a probability 
>90% of lying at 2 > 8 after the inclusion of the flux 
constraint (0±12 nJy: nominally the flux constraint one 
would obtain for a z ~ 9-10 galaxy based on a single 
orbit of Yios-band observations) and have a measured 
Tl 25 —-ffi 60 color >0.5 mag. Our simulation results make 
it clear how important the pre-selection can be. While 
increasing the efficiency of our search results significantly, 
pre-selection can also introduce a modest amount of in¬ 
completeness into the z ~ 9-10 samples we identify from 
CANDELS, particularly at z < 9 (where it is ~40% from 
the pre-selection step alone). 

In addition to our considering the selection of 
sources from CANDELS-UDS, CANDELS-COSMOS, 
and CANDELS-UDS fields, we also consider the selec¬ 
tion of z ~ 9 and z ~ 10 galaxies from the CAN- 
DELS GOODS-North, CANDELS GOODS-South, and 
ERS fields. 

In computing the number of confirmed z = 9-10 
sources from our program, we assume all the sources in 
Table @] are bona-fide z ~ 9-10 galaxies and there is no 
contamination in our selection. This would appear to be 
a good assumption, given that the typical z ~ 9-10 candi¬ 
date formally prefers a z > 8 solution at 99% likelihood. 
We do not include EGS910-2 and UDS910-1 in our LF 
calculation, since they do not meet our formal criteria for 
inclusion (but nevertheless appear to be probable z > 8.5 
galaxies). We suppose that all of the candidates from our 
follow-up program that were not explicitly confirmed by 
that program to lie at z < 8.4 (all but one of these can¬ 
didates was detected at >2er in the follow-up Yios-band 
observations and are therefore unlikely z > 8.4). The 
0 ~ 9 candidate GS-z9-5 is modeled as 0.5 z ~ 9 galaxy 
(i.e., ~50% probability of contamination), given its com¬ 
puted P(z > 8) was only 0.66 (Figure [SJ. We ignore the 
impact of possible lensing amplification on one source in 
our selection (EGS910-3) given the size of the magnifica¬ 
tion factor (0.25 mag) and the fact that source volume 
and magnification factor trade off in such a way to have 
little impact on the derived LF. 

Our z ~ 9 and z ~ 10 LF results are presented in 
Figure [7] and also tabulated in Table [5] In computing 
the uncertainties on the z ~ 9 and z ~ 10 UV LFs, we 
also include the expected large-scale structure uncertain¬ 
ties, using the results from the cosmic variance calculator 
of Trenti & Stiavelli (2008) and the observed comoving 
volume density. For context, the earlier LF results of 



Fig. 9.— The present determinations of the UV luminosity den¬ 
sity of galaxies at z ~ 9 and z ~ 10 ( solid red circles ) brightward of 
—20 mag using the present search over all 5 CANDELS fields. The 
black crosses show the luminosity density of galaxies brightward of 
—20 mag determined by Bouwens et al. (2015) at z ~ 4, z ~ 5, 
z ~ 6, z ~ 7, and z ~ 8. The shaded region shows the evolution¬ 
ary trend to z > 8 in the UV luminosity density preferred at 68% 
confidence. The red line shows the observed evolution in the UV 
luminosity density beyond z ~ 8 and can be approximately repre¬ 
sented by dlog 10 puv/dz = —0.45 d= 0.07 dex. The UV luminosity 
density we find for luminous z ~ 9 and z ~ 10 galaxies lies below 
a simple extrapolation of the evolution from z ~ 4 to z ~ 8. 

McLure et al. (2013), Oesch et al. (2013, 2014), Bouwens 
et al. (2014b, 2015), and McLeod et al. (2015) are pre¬ 
sented in Figure [7] 

Given the limited numbers of z ~ 9-10 candidates in 
our samples and some arbitrariness in the choice of bin 
centers (and width) for our stepwise 2 ~ 9-10 LF results, 
it is conceivable that our z ~ 9-10 LF results could de¬ 
pend on how we bin our sample. We therefore also model 
the bright end of the LF as a power law (motivated by 
the results of e.g. Bowler et al. 2014 and Bouwens et al. 
2015). By marginalizing over both the normalization and 
power-law slope of the model LFs, we derive constraints 
( 68 % confidence levels) on the volume density of z ~ 9- 
10 galaxy candidate for a given Mjjv ■ These results are 
presented in Figure [7] and shown in Table [ 6 ] 

To guide expectations, we also included on Figure [7] 
the LFs we derive extrapolating the LF results from 
Bouwens et al. (2015) to * ~ 9 and z ~ 10.2 (the mean 
redshift of our z ~ 9 and z ~ 10 derived from our se¬ 
lection volume simulations). In deriving a fitting for¬ 
mula from the Bouwens et al. (2015) results, we only 
consider the evolution over the range z ~ 8 to z ~ 5 
where the evolution of the UV LF evolves in a rela¬ 
tively smooth manner. This results in the Schechter pa¬ 
rameters depending on redshift in the following manner: 
Miry = (-20.97 ± 0.10) + (0.17 ± 0.10)(z - 6 ), <j>* = 
(0.45t[);^)10(- o ' 21 ±o - o 9 ^ 2 - 6 h0- 3 Mpc _3 , a = (-1.91 ± 
0.05) + (—0.13±0.05)(.z — 6 ). This fitting formula implies 
a M* = -20.45, <j>* = 0.10 x 10 “ 3 Mpc“ 3 , and a = -2.3 
at 2 ~ 9 and M* = -20.28, 0* = 0.059 x 10 “ 3 Mpc" 3 , 
and a = —2.46 at z ~ 10.2. The best-fit evolutionary sce¬ 
nario from Bowler et al. (2014) is very similar to what we 
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present above (where it was found that cLM*/dz ~ 0 . 2 ). 

Relative to these extrapolations of z = 4-8 results to 
z > 8 , our present LFs are typically 1.5 x lower in the 
mean, consistent with a slightly faster evolution at z > 8 . 
Results from the above fitting formula are also featured 
on Figure [ 8 ] where we combine the new LF results with 
previous results from the literature at 2 ~ 4, z ~ 5, 
z ~ 6, z ~ 7, and z ~ 8 . 

The present estimates of the volume density of par¬ 
ticularly bright z ~ 9 galaxies are the first available in 
the literature and have comparable uncertainties to the 
bright end of the present LF determinations at z ~ 10 
(and as derived earlier: Oesch et al. 2014; Bouwens et al. 
2015). 

7.2. Evolution of the UV Luminosity Density for 
Luminous Galaxies from z ~ 10 to z ~ 4 

With the present constraints on the volume density of 
bright galaxies at 2 ~ 9 and z ~ 10, we can now examine 
the evolution of the luminosity density of galaxies in the 
rest-frame UV with cosmic time. 

Integrating up the light in our binned representation of 
the bright end of the z ~ 9 and 2 ~ 10 LFs to —20 mag 
(approximately the faint-end limit of the present bright 
search), we derive a total UV luminosity density puv 

of IO 24 ' 21 - 019 erg/s/Hz/Mpc 3 at 2 ~ 9 and io 2394 -o .28 
erg/s/Hz/Mpc 3 at 2 ~ 10. These inferred luminosity 
densities are 5^ x and 8+3 x lower, respectively, than 
Bouwens et al. (2015) found at 2 ~ 8 . 

In Figure [9l we compare the derived luminosity density 
puv at z ~ 9 and 2 ~ 10 with the UV luminosity density 
Bouwens et al. (2015) derived from 2 ~ 4 to 2 ~ 8 to the 
same limiting magnitude. Also shown on this figure is the 
best-fit evolution ( 68 % confidence intervals) derived for 
the UV luminosity density based on the 2 ~ 4-8 results. 
Our new luminosity density results puv at 2 ~ 9 and 2 ~ 
10 are 2 . 6 + 5 'g x and 2 . 2+ 2 ' 9 x lower than the extrapolated 
luminosity densities at 2 ~ 9.0 and 2 ~ 10.2. 

The present results are therefore consistent with a 
more rapid evolution at 2 > 8 than between 2 ~ 8 
and 2 ~ 4, as first noted by Oesch et al. (2012), with 
a best-fit d log 10 puv/dz evolution of —0.45 ± 0.07 vs. 
the —0.29 ± 0.02 evolution observed from 2 ~ 8 to 2 ~ 4. 
Interestingly enough, the evolution we derive here is com¬ 
pletely consistent with the d log 10 puv/dz = — 0.54l o ;gg 
scaling that Oesch et al. (2012) previously derived. It is 
also consistent with the (1 + z)~ 10 9 evolutionary scal¬ 
ings considered in Oesch et al. (2014: see also Oesch et 
al. 2013). 

Other authors (Ishigaki et al. 2015; Oesch et al. 2015; 
McLeod et al. 2015) have remarked that there is less ev¬ 
idence for a faster evolution in the UV luminosity den¬ 
sity at 2 > 8 integrated to lower luminosities, particu¬ 
larly when incorporating new results from the full Hub¬ 
ble Frontier Fields program (Lotz et al. 2014; Coe et al. 
2015). If the HUDF/XDF is systematically underdense 
in z ~ 9-10 galaxies or if the evolution of the LF pro¬ 
ceeds more smoothly for lower luminosity galaxies, one 
might expect such a result. A much more definitive ex¬ 
ploration of this issue will be possible by considering the 
much larger number of faint sources at 2 ~ 9-10 expected 
from the full Frontier Fields program (Coe et al. 2015). 


8 . SUMMARY 

In this paper, we provide improved constraints on the 
volume density of especially luminous [H\qq,ab < 27) 
2 ~ 9 galaxies by making use of a search over all 5 CAN- 
DELS fields (750 arernin 2 area in total). We also rederive 
the bright end of the 2 ~ 10 LF and extend it slightly 
fainter taking advantage of the additional search power 
provided by the S-CANDELS data (Ashby et al. 2015) 
over the CANDELS-UDS, COSMOS, and EGS fields. 

To obtain these constraints, we extend the earlier sam¬ 
ple of 6 bright 2 ~ 9 and 2 ~ 10 galaxies identified 
by Oesch et al. (2015) to make use of the search area 
over all 5 CANDELS fields, including an additional ~ 
450 aremin 2 area over the CANDELS-UDS, CANDELS- 
COSMOS, and CANDELS-EGS fields. We also expand 
the Oesch et al. (2014) search over the GOODS-South 
and GOODS-North to obtain a more complete selection 
of galaxies over the redshift range 2 ~ 8 .4-9.0 (see §5). 

To identify a robust sample of bright 2 ~ 9 and 2 ~ 10 
galaxies over this new area, we employed a two-part 
strategy for selecting sources. To begin, we selected 
those sources that showed the highest probability of cor¬ 
responding to 2 ~ 9 and 2 ~ 10 galaxies based on the 
existing observations. Such identifications were possible, 
given the depth of the HST and ground-based imaging 
observations which placed constraints on both the sharp¬ 
ness of a possible spectral break at 1.2pm and also the 
spectral slope blueward of the break. 

The second step was to obtain deep observations on 
each of these candidates at 1.05/zm to test the nature 
of these candidates. These follow-up observations were 
obtained with the 11-orbit HST program z9-CANDELS 
(Bouwens 2014: GO 13792). 

Using the new HST observations from the z9- 
CANDELS program, we find that 5 out of the 12 2 ~ 9-10 
candidates we have followed up with our program appear 
to be bona-fide bright 2 ~ 9-10 galaxies (Figure H}. 

Combining our new samples with previous samples of 
luminous 2 ~ 9-10 candidates over CANDELS GOODS- 
North and GOODS-South (Oesch et al. 2014) and also 
including four probable 2 ~ 8 .4-9.0 galaxies we identified 
over the GOODS-S field (finding no additional bright 
sources over the GOODS-N field), we identify a total 
sample of 15 bright 2 ~ 9-10 candidates over the five 
CANDELS fields (Table [4}. This is >2 x larger than the 
sample of luminous 2 ~ 9-10 galaxies we had identified 
earlier in Oesch et al. (2014). 

We use these larger samples of 2 ~ 9-10 galaxies to de¬ 
rive improved constraints on the bright end of the 2 ~ 9 
and 2 ~ 10 LFs. We achieve these results by simulating 
both stages in the selection process (for the CANDELS- 
UDS, COSMOS, and EGS fields) - or a single stage in the 
case of the CANDELS GOODS-North, GOODS-South, 
or ERS fields - to arrive at statistically-robust conclu¬ 
sions. 

As one would expect with any follow-up program of 
limited size, the z9-CANDELS program is not able to 
observe all potential 2 ~ 9 and 2 ~ 10 galaxies over the 
CANDELS-UDS, COSMOS, and EGS fieldsQand there- 
fore likely suffers from moderate incompleteness, partic- 

12 Several possible examples of missed z ~ 9-10 candidates could 
include those lower-probability candidates tabulated in Appendix 
C. 
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ularly for galaxies at z < 9, due to our exclusively pre¬ 
selecting for follow-up sources with redder J 125 — i?i 6 o 
>0.5 colors than galaxies which exist in this range. Some 
incompleteness would also result from the detection of 
spurious flux in the optical bands (~25% effect). Nev¬ 
ertheless, we remark that all such effects are included 
in our selection volume simulations and our focusing on 
J 125 — > 0.5 galaxies means our search is most com¬ 

plete at z > 9. 

We would expect these results to be improved, particu¬ 
larly towards fainter magnitudes, if even deeper observa¬ 
tions over the CANDELS-UDS, CANDELS-COSMOS, 
and CANDELS-EGS fields could be obtained. Deeper 
observations in the F105W band should increase the size 
of our bright (-Hi 6 o,ab < 27) z ~ 9-10 samples over the 
CANDELS-WIDE fields by improving the completeness 


of our search results at z < 9(3 Meanwhile, deeper ob¬ 
servations in various redder bands (e.g., JHuq^ab) would 
allow us to extend these searches to even fainter magni¬ 
tudes (i.e., to >26.5 mag). We would also expect gains 
from analysis of the 500-orbit, 500 arcmin 2 BoRG[ 2 g 10 ] 
program (Trenti 2014). 

We thank Jim Dunlop and Steve Finkelstein for valu¬ 
able conversations. We thank Steve Willner for pro¬ 
viding us with feedback on an earlier draft of this 
manuscript. This work is based in part on observations 
made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. We acknowl¬ 
edge the support of NASA grant NAG5-7697, NASA 
grant HST- GO-11563, and a NWO vrij competitie grant 
600.065.140.11N211. 
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APPENDIX 

A. POSTAGE STAMPS AND SEDS OF THE TARGETED CANDIDATES 

We identified 12 candidate z = 9-10 sources over the CANDELS-UDS, CANDELS-COSMOS, and CANDELS-EGS 
fields that we explicitly targeted for follow-up observations. 

Postage stamp images of these 12 candidates are presented in Figure [10] As should be obvious from this figure, all 
six of the present z ~ 9-10 candidates show clear detections in the i?i 6 o-band, as well as significant ~ 2-3 ct detections 
in the Ji 25 -band and J-ffi 4 o-band observations (where available), as well as in the S-CANDELS Spitzer/IRAC data. 

The observed photometry to those candidate sources are presented in Figure [TT| along with SED fits to a model 
z > 8 galaxy and a model z < 3 galaxy. Also shown on this figure is the redshift likelihood distribution ( solid black 
line) based on the photometry we have available for each candidate in the ~20 different wavelength channels ( HST + 
Spitzer/IRAC + ground-based observations). In addition, this figure presents the redshift likelihood distribution we 
would expect, assuming these candidates are not detected in the single orbit of follow-up Yios-band observations from 
the z9-CANDELS program. 

All 12 candidates show similar Hi 60 — [3.6] colors (see also Wilkins et al. 2016), red [3.6] — [4.5] colors, and observed 
sizes (Holwerda et al. 2015; Shibuya et al. 2015) to the first samples of particularly luminous 2 ~ 9-10 galaxies identified 
by Oesch et al. (2014). 

B. TARGETED Z ~ 9-10 CANDIDATES THAT ARE MOST LIKELY AT Z < 8.4 

A fraction of the candidate z ~ 9-10 galaxies targeted by our follow-up observations from the z9-CANDELS program 
appear not to be at z > 8 . In Figure [l2l we present postage stamps for 9 such sources (7 explicitly targeted by our 
program and 2 lower probability z ~ 9-10 candidates which were incidentally targeted). In almost all of the candidates 
which are not confirmed by our follow-up observations, the sources show a > 2 er detection in the Yios-band. Figure [13] 
shows the best-fit SED models for the tentative z ~ 9-10 galaxies that were not confirmed by observations from our 
follow-up program. 

Detailed comments on specific candidate z ~ 9-10 galaxy that were targeted for follow-up with our z9-CANDELS 
program can be found below: 

COS910-0: Follow-up of COS910-0 in the Yios-band shows a 3er detection at 1.05/im. A detailed fit to its SED 
suggests that it is actually a star-forming galaxy at 2 = 7.8. Its inclusion in our original sample of z > 8.4 candidate 
galaxies occurred due to its measured J 125 — -H 16 O color, which was likely redder than reality due to the impact of 
noise. 

COS910-2: COS910-2 is detected at 2.2er in the Yios-band follow-up observations. Such a detection is not expected 
if the galaxy is actually at z > 9, and so the source must have a redshift of z < 8 . It may correspond to either a 
0 ~ 0.5 or a z ~ 8 galaxy. 

COS910-3: Follow-up of COS910-2 in the Yios-band and the JiJi 4 o-band (1/3 of an orbit) shows a la detection at 
1.05/im and 2 a detection in the Ji 7 i 4 o band. Perhaps, most importantly, the deeper photometry over the source in 
the JH 140 band confirms that the source has a Ji?i 4 o — Hi 60 color of 1.25 mag. While this is consistent with the 
source having a redshift of the z > 11 , one would not expect to detect the source at 2 er in the J 125 band in this case. 
These results suggest that the apparent break in the spectrum at 1.5/zm is not especially sharp and this source shows 
faint but detectable flux to much bluer wavelengths. 

UDS910-0: Follow-up of UDS910-0 in the Y 105 band shows a 4<r detection at 1.05/im. Therefore, this source cannot 
plausibly correspond to a z > 8.4 galaxy. The best-fit redshift we compute for the source is z ~ 7. Like COS910-0, its 
inclusion in our sample of 2 ~ 9-10 galaxies likely occurred as a result of noise in the measured J 125 — Hwo color. 
EGS910-1: Follow-up of EGS910-1 in the Yios-band shows a 3er detection at 1.05/rm, consistent with a redshift of 
0 < 8. 
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Fig. 10. — HST + Spitzer /IRAC images (6 ;/ x 6") of all preselected candidate z ~ 9-10 galaxies over the CANDELS-UDS + CANDELS- 
COSMOS + CANDELS-EGS fields that were preselected for targeted follow-up with the z9-CANDELS program. The best-fit model flux 
from neighboring sources has been removed in the Spitzer /IRAC images shown here. Our preselected CANDELS-UDS and CANDELS- 
COSMOS candidates (1) are estimated to show a >50% probability of corresponding to a z > 8 galaxy and (2) can be confirmed to lie at 
z > 8 with >90% confidence with the addition of a single orbit of HST follow-up observations (assuming a flat redshift prior and the EAZY 
SED template set). Each of the preselected CANDELS-EGS candidates shown here can be confirmed to be secure at >90% probability with 
a single orbit of HST follow-up observations. Each of these candidates was subject to 1-orbit follow-up observations with our z9-CANDELS 
program at 1.05/xm. 
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Fig. 11.— (leftmost column ) Best-fit SED models to the observed RST+Spitzer/IRAC-l-ground-based photometry of the preselected 
candidate z ~ 9-10 galaxies we have identified for targeted follow-up with the z9-CANDELS follow-up program. The dotted red upper 
limits show the approximate constraints we will be able to set on the 1.05//rn fluxes of the candidates, assuming they are at z > 8.4. The 
points and lines are otherwise as in Figurc[4] (second leftmost column) Redshift likelihood distribution for the same 12 z ~ 9-10 candidates 
using current observations (solid lines) and also making use of our 1 -orbit Y) 05 -band follow-up observations (dotted lines ) assuming the 
sources are at 2 > 9. P(z < 8 ) and P new {z < 8 ) indicate the probability that our candidates are at 2 < 8 using the current observations 
and including our follow-up observations (again assuming they are 2 > 9). It should be clear that our z9-CANDELS follow-up observations 
should significantly improve our confidence in the present set of 2 > 9 candidates, increasing it from 72-95% to 95.8-99.8%. 


EGS910-4: Follow-up of EGS910-4 in the Yios-band shows a 3er detection at 1.05/tm, consistent with a redshift of 
z < 8. 

EGS910-5: Follow-up of EGS910-5 in the y 10 5 -band shows a 1.7er detection at 1.05/itm, strongly suggesting it is not 
at z > 8. 

EGS910-6: EGS910-6 is detected at 4.5er in our y 10 5 -band follow-up observations at 1.05/um, providing clear evidence 
it is at z < 8. A 2 < 8 solution is preferred at 99.8% confidence iFigurefl^l). 

EGS910-7: Follow-up of EGS910-7 in the Yios-band only shows a l.ler detection at 1.05/rm. However, the Hiqq — [3.6] 
color is sufficiently red that it seems more consistent with a lower-redshift galaxy than az>8 galaxy. The nature of 
this source is unclear. 
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Fig. 12.— HST + Spitzer/IKAC images for 7 targeted candidate z ~ 9-10 galaxies which follow-up observations from our HST z9- 
CANDELS program failed to confirm as probable z > 9 galaxies using HST follow-up observations with our z9-CANDELS program. 
F105W-band observations were also obtained for the sources shown in the lowest 2 rows of this figure from the z9-CANDELS HST 
program. While neither source was not explicitly targeted for observations by our program due to a relatively low prior probability of being 
at z > 8, we could not rule out that possibility and so took advantage of the new data to gain more insight into their nature. 















22 



log l0 Wavelength (nm) 


Redshif t 


Fig. 13.— {left) Best-fit SED models to the observed ifS'T+S'piteer/IRAC+ground-based photometry for 7 targeted candidate z ~ 9-10 
galaxies (COS910-0, COS910-2, COS910-3, EGS910-1, EGS910-4, EGS910-5) which were not confirmed as probable z > 8.4 galaxies using 
follow-up observations from the z9-CANDELS program. The points and lines are otherwise as in Figure[4] F105W-band observations were 
also obtained for the sources shown in the lowest 2 rows of this figure (EGS910-6 and EGS910-7). While neither source was not explicitly 
targeted for observations by our program due to a relatively low prior probability of being at z > 8, we could not rule out that possibility 
and so took advantage of the new data to gain more insight into their nature, [right) Redshift likelihood distribution for these z ~ 9-10 
candidates incorporating both our follow-up observations with HST and the TfST+iSpitzer/IRAC+ground-based observations that were 
used in the pre-selection (solid lines). 


C. OTHER CANDIDATE Z ~ 9-10 GALAXIES OVER THE CANDELS-UDS, COSMOS, AND EGS FIELDS 

When putting together a compilation of the most promising z ~ 9-10 candidates to target with follow-up HST 
observations, we experimented with a variety of different procedures to identify possible z ~ 9-10 galaxy candidates. 
As a result, we identified a larger number of possible z = 9-10 galaxy candidates than we could thoroughly follow up 
with the 11 orbits allocated to our HST program. 

While the redshift likelihood distributions we derived for these candidates suggested that the vast majority of them 
likely corresponded to sources at 2 < 8, it is nevertheless possible that several of these candidates might have redshifts 
in excess of z ~ 8. 

We provide a short list of the “lower-probability” z ~ 9-10 galaxy candidates we identified over the CANDELS 
COSMOS, UDS, and EGS fields in Table [7] These sources did not meet our criteria for preselection, but nevertheless 
may correspond to z ~ 9-10 galaxies. Figure [14] shows postage stamp cut-outs of these sources, while Figure [15] 
illustrates fits to their spectral energy distributions and our derived redshift likelihood contours for these sources. 
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TABLE 7 

Possible z ~ 9-10 Galaxies over the CANDELS UDS, COSMOS, and EGS program that did not satisfy for our criteria for 

PRESELECTION AND HENCE WERE NOT TARGETED FOR FOLLOW-UP OBSERVATIONS. 


ID 

R.A. 

Dec 

# 160 , AB 

Zphot 

V 

00 

COS910-4 

10:00:15.52 

02:17:01.5 

25.6±0.1 

9.3 

0.18 

UDS910-2 b 

02:17:13.08 

-05:15:55.4 

26.6±0.2 

10.2 

0.68 

UDS910-3 

02:17:52.38 

-05:15:06.3 

26.9±0.2 

9.4 

0.28 

UDS910-4 

02:17:14.61 

-05:15:15.7 

26.6±0.2 

9.1 

0.50 


a Redshift likelihood is computed using the flux measurements for these sources in all photometric bands shown in Table [T] 
b While the source UDS910-2 nominally has >50% probability of lying at z > 8, it was not targeted with our follow-up program z9- 
CANDELS because it could not be confirmed to be >90% probability candidate with the addition of a single orbit of HST observations. 
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Fig. 14.— HST + Spitzer /IRAC images (6 ;/ x O 77 ) of 4 possible z ~ 9-10 candidate galaxies that did not meet our crit eria for preselection 
and hence were not targeted by our z9-CANDELS program. Sources are shown in the same passbands, as in Figure fTol None of these 
sources is explicitly targeted with our z9-CANDELS follow-up program. 
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Fig. 15.— (left) Best-fit SED models to the observed HS T+ Spitzer /IRAC+ground-based photometry of 4 possible 2 ~ 9-10 candidate 
galxies that did n ot m eet our criteria for preselection and hence were not targeted by our z9-CANDELS follow-up program. Symbols are 
as in Figure [4] and fTTI (right) Redshift likelihood distribution for these 6 9-10 candidates. 





































































